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The lattice cluster theory (LCT) for semiflexible linear telechelic melts, developed in paper I, is applied 
to examine the influence of chain stiffness on the average degree of self-assembly and the basic thermody¬ 
namic properties of linear telechelic polymer melts. Our calculations imply that chain stiffness promotes 
self-assembly of linear telechelic polymer melts that assemble on cooling when either polymer volume fraction 
(j) or temperature T is high, but opposes self-assembly when both </> and T are sufficiently low. This allows us 
to identify a boundary line in the (fyT plane that separates two regions of qualitatively different influence of 
chain stiffness on self-assembly. The enthalpy and entropy of self-assembly are usually treated as adjustable 
parameters in classical Flory-Huggins type theories for the equilibrium self-assembly of polymers, but they 
are demonstrated here to strongly depend on chain stiffness. Moreover, illustrative calculations for the depen¬ 
dence of the entropy density of linear telechelic polymer melts on chain stiffness demonstrate the importance 
of including semiflexibility within the LCT when exploring the nature of glass formation in models of linear 
telechelic polymer melts. 


I. INTRODUCTION 

Despite the increasing scientific interest and techno¬ 
logical importance of telechelic polymers, a molecular 
theory that predicts the relation between the monomer 
molecular structure and the thermodynamic properties 
in such systems has been slow to develop. Prior the- 
oriesi^— of self-assembly in telechelic polymers tradi¬ 
tionally employ highly coarse grained models that rep¬ 
resent the assembling molecular species as a structure¬ 
less entity and hence, cannot address the influence of lo¬ 
cal molecular structure on the self-assembly of telechelic 
polymers. Dudowicz and Freedii have recently reformu¬ 
lated the lattice cluster theory (LCT) for the thermody¬ 
namics of polymer systems to treat strongly interacting, 
self-assembling polymers composed of fully flexible linear 
telechelic chains. Because the LCT employs an interme¬ 
diate level of coarse grained models that retain minimal 
aspects of monomer molecular structure and interactions 
in telechelic polymers, the theory enables establishing the 
relation between the molecular structure dependent in¬ 
teraction parameters of the model and the thermody¬ 
namic properties of these complex fluids i^^d^ Therefore, 
the LCT provides a potential framework for obtaining in¬ 
formation that is useful for designing telechelic polymer 
materials. 

Paper I of this series^^ further extends the LCT for the 
thermodynamics of linear telechelic polymers to include 
a description of chain semiflexibility. The significance of 
this extension is at least two-fold. Firstly, computer sim- 
ulationsi^^— indicate that chain semiflexibility strongly 
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affects the self-assembly process and structural proper¬ 
ties of telechelic polymers. The LCT for semiflexible 
telechelic polymers provides a theoretical tool for inves¬ 
tigating the influence of chain stiffness on the thermody¬ 
namic properties of self-assembling telechelic polymers 
and hence, offers the possibility of better understand¬ 
ing experiments and simulations of telechelic polymers. 
Secondly, chain semiflexibility is crucial for exploring the 
behavior of glass formation in polymers within the gen¬ 
eralized entropy theory (GET)ji^ a theory of polymer 
glass formation that merges the LCT for the thermody¬ 
namics of semiflexible polymers^^ with the Adam-Gibbs 
(AG) relation between the structural relaxation time and 
the configurational entropyr2i2i because models of fully 
flexible polymer chains fail to exhibit the characteristic 
glassy behavior in the GET, as noted in Ref. [l^. 

The previous wor k^^d^ analyzes the influence of vari¬ 
ous molecular parameters on the thermodynamic prop¬ 
erties only for fully flexible linear telechelic polymers. 
The influence of chain semiflexibility upon the thermo¬ 
dynamics of self-assembling telechelic polymers is inves¬ 
tigated using the theoretical advances of the LGT in pa¬ 
per to explore the relation between molecular fac¬ 
tors and thermodynamic properties in equilibrium self¬ 
assembling semiflexible telechelic polymers, and in par¬ 
ticular, in semiflexible linear telechelic polymers. 

Section II provides a brief summary of the LCT for 
semiflexible linear telechelic polymers. Section III be¬ 
gins by examining the influence of chain stiffness on the 
average degree of self-assembly of linear telechelic poly¬ 
mers that assemble on cooling. Our calculations demon¬ 
strate that chain stiffness can either promote or oppose 
the self-assembly of linear telechelic polymer melts, de¬ 
pending on the polymer volume fraction and the temper¬ 
ature. Specifically, while chain stiffness promotes self- 
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FIG. 1. Schematic illustration of the lattice model for a self- 
assembled linear telechelic cluster composed of three short 
chains with M = 5 united atom groups in each chain. Solid 
circles designate the chain’s ends that can participate in sticky 
interactions, while open circles denote the united atom groups 
in the chain interior. Solid lines represent the chemical bonds 
between two united atom groups within the same chain, while 
dashed lines denote the sticky physical bonds between two 
stickers. The model prescribes two different nearest neighbor 
interaction energies e and ts for ordinary and sticky-sticky 
interactions, respectively. 
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assembly of linear telechelic polymer melts when either 
polymer volume fraction (jj or temperature T is high, self- 
assembly is suppressed by chain stiffness when both ip 
and T are sufficiently low. This allows us to identify a 
boundary line in the (p-T plane, which separates two re¬ 
gions in which the influence of chain stiffness on the self- 
assembly is qualitatively different. We further show that 
the boundary line is independent of the sticky interac¬ 
tion energy but changes with bending energy, molecular 
weight, and van der Waals interaction energy. Section III 
then presents calculations of the enthalpy and entropy of 
self-assembly that are usually treated as adjustable pa¬ 
rameters in classical Flory-Huggins (FH) type theories 
for the equilibrium self-assembly of polymers (e.g., see 
Refs. [23I - I2CT but that are shown here to depend strongly 
on chain stiffness. Section III further presents calcula¬ 
tions for the influence of chain stiffness on the entropy 
density of linear telechelic polymer melts, thereby illus¬ 
trating the importance of including semiflexibility within 
the LCT for exploring the nature of glass formation in 
models of linear telechelic polymer melts. 


II. LATTICE CLUSTER THEORY FOR SEMIFLEXIBLE 
LINEAR TELECHELIC POLYMER MELTS 

We briefly review the basic characteristics of the LCT 
for semiflexible linear telechelic polymer melts to pro¬ 
vide necessary definitions, notation, and description of 
the model; more details can be found in paper L— 

The system consists of m semiflexible linear chains that 
are located on a simple cubic lattice with lattice coordi¬ 
nation number equal to z = 6 and with Ni total lattice 
sites. Since the system is treated as a compressible melt, 
the lattice also contains empty sites whose volume frac¬ 
tion is determined from the equation of state for a system 
at constant temperature and pressure. The volume frac¬ 
tion of polymer chains is defined hy (p = mM/Ni, where 
M is the number of united atom groups per chain. In 
order to model telechelic chains, each chain’s end seg¬ 


ments (called “stickers” and represented as solid circles 
in Fig. 1) are distinguished from the other united atom 
groups lying in the chain interior (depicted by open cir¬ 
cles in Fig. I). Two stickers are allowed to form a sticky 
“bond” (pictured as dashed lines in Fig. I) and interact 
with an attractive sticky interaction energy Cg when they 
are located on nearest neighbor lattice sitesjii A nega¬ 
tive Cs implies that polymeric clusters form upon cooling. 
This sticky energy eg may greatly exceed the microscopic 
van der Waals interaction energy e, a parameter that de¬ 
scribes the attractive interaction strength between two 
non-stickers as well as between a sticker and a non-sticker 
(see Fig. 1) and that is conventionally treated as posi¬ 
tive in the original LCTi21i2^ The current LCT assumes 
that the stickers are mono-functional (i.e., the stickers at 
each of the ends of the telechelics can only participate 
in one sticky interaction) and that both cyclic and linear 
associative clusters may form, in accord with the pre¬ 
vious worbiii^^ and the analysis of Jacobson and Stock- 
mayei'i^ Following the previous treatment in the LCT^ 
chain semiflexibility is treated by introducing a bending 
energy penalty Eh whenever a pair of consecutive bonds 
from a single telechelic chain lies along orthogonal direc¬ 
tions. Because the sticky bonds are treated as being fully 
flexible, the theory leads to a model where the associ¬ 
ated chains are monomeric rods with fully freely rotating 
junctions at the sticky bonds in the limit Eh ^ oo. This 
feature prevents the occurrence of crystallization in rigid 
linear telechelic polymers, and hence, may be important 
for exploring glass formation in such systems. 

As derived in paper I^ the Helmholtz free energy / 
per lattice site of a semiflexible telechelic polymer melt 
is the sum of the free energy fo of the hypothetical ref¬ 
erence system in the absence of sticky interactions and 
the free energy contribution fg arising from the sticky 
interactions. 


l3f = P3fo + Pfs, ( 1 ) 

where /3 = l/{kBT) with ks being Boltzmann’s constant 
and T designating the absolute temperature. The free 
energy fo of the hypothetical reference system is identical 
to the system under consideration, except that the sticky 
interaction is absent. Thus, the free energy fo of the 
reference system is independent of Cg and its derivation 
is described elsewhere in detai h^^i^° 

The sticky contribution fg is derived by hrst defining 
the density y of sticky bonds (i.e., the ratio of the num¬ 
ber of sticky bonds in the system to the total number of 
lattice sites) and emerges as 

4 
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where x = 2/M is the fraction of stickers in a single chain, 
and Yi {i = are corrections to the zeroth-order 

mean-field free energy arising from short range 

correlations up to the scale of four consecutive bonds. 
Explicit expressions for (i = are provided in 

paper The variable y in Eqs. (2) and (3) is deter¬ 
mined by the maximum term method?^ i.e., by applying 
the condition, 


dWs) 

dy 


= 0 . 

Ni,t,4> 


( 4 ) 


The solution of Eq. (4), i.e., the concentration y* of 
the sticky bonds, is then substituted into Eqs. (1-3), 
leading to the final expression for the free energy / of a 
semiflexible linear telechelic melt, 


13f =13fo - (l)x\n{(j)x) + {(j)x - 2y*) ln((/)a; - 2y*) 
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-J2Y{y*y, (5) 


which depends on all molecular and thermodynamic pa¬ 
rameters including T, (j), M, e, Ef,, and e^. Notice that 
the volume fraction of the active stickers (i.e., those par¬ 
ticipating in sticky interactions) is simply 2y* and that 
the upper limit for y* is = <^/M, because the 

chain’s ends are each assumed to be mono-functional in 
the present model. 

While the present paper discusses the model and re¬ 
sults with reference to compressible melts, the mathe¬ 
matical equivalence between the excess thermodynamic 
properties of a compressible melt and an incompressible 
solution enables describing both types of systems. The 
model for an incompressible solution consists of one-bead 
solvent molecules replacing the empty lattice sitesjii The 
free energy of a compressible polymer melt is isomor¬ 
phic to that of an incompressible polymer solution, with 
the e parameter being replaced by the exchange energy 
Eex = Cpp + ess — 2eps, where e^p, Css and Cps represent the 
strengths of the nearest neighbor interaction between two 
polymer segments, two solvent molecules and a polymer 
segment and a solvent molecule, respectively. 


III. RESULTS AND DISCUSSION 


This section examines the influence of chain stiffness 
on the average degree of self-assembly and the enthalpy 
and entropy of self-assembly of linear telechelic polymer 
melts. This section also demonstrates the importance 
of including chain semiflexibility within the LCT for ex¬ 
ploring glass formation in the model of linear telechelic 
polymer melts. 



FIG. 2. i/* as a function of bending energy Eb for var¬ 

ious sticky interaction energies es- The computations are 
performed for a melt of linear telechelic chains, where the 
polymer volume fraction is (f) = 0.9, the molecular weight is 
M = 100, and the van der Waals interaction energy is e = 100 
K. The temperature is fixed at T = 200 K. 
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FIG. 3. y* as a function of bending energy Eb for various 
molecular weights M. The computations are performed for 
a melt of linear telechelic chains, where the polymer volume 
fraction is 0 = 0.9, the van der Waals interaction energy is 
e = 200 K, and the sticky interaction energy is Sa = —100 K. 
The temperature is fixed at T — 200 K. 


A. Influence of chain stiffness on the average degree of 
self-assembly 


An understanding of the self-assembly of telechelic sys¬ 
tems is crucial for many practical applications^ir— As 
explained in Sec. II and paper I^ the concentration 
y* of the sticky bonds measures the volume fraction of 
the active stickers in the LCT for linear telechelic poly¬ 
mer melts. Hence, y* directly reflects the average degree 
of self-assembly and other important quantities such as 
the order parameter <I> and the transition temperature 
Tp of self-assembly3 (The average degree < A^ > of self- 
assembly is more precisely defined by < N 1/(1 — 4>) 
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FIG. 4. y* as a. function of bending energy Eb for various poly¬ 
mer volume fractions (j> for a melt comprising linear telechelic 
chains, where the molecular weight is M — 100, the van der 
Waals interaction energy is e = 200 K, the sticky interac¬ 
tion energy is £s = —100 K, and the temperature is fixed at 
T = 200 K. 



FIG. 5. Contour plot of y* in the (f>-T plane. The compu¬ 
tations are performed for a melt of linear telechelic chains, 
where the molecular weight is M = 100, the van der Waals 
interaction energy is e = 200 K, the sticky interaction energy 
is £s = —100 K, and the bending energy is Eb = 2000 K. 
The solid line denotes the boundary marking the states with 
J/‘ =0. 

with $ = y*/y^ax'^ see Ref. [l^ for more details.) We thus 
focus on the general trends exhibited by the variation of 
y* with bending energy. Because the dependence of y* on 
various molecular and thermodynamic parameters, such 
as M, e, and e^, has been extensively examined in 

Refs. [l^lT^ for the model of fully flexible linear telechelic 
polymers, the computations presented in this subsection 
are devoted to understanding the role of chain stiffness in 
the self-assembly of linear telechelic polymer melts. As 
in previous work j^^d^ the lattice coordination number is 
taken as z = 6 in all our computations. 

Because the quantitative effect of Eb on y* is generally 


found to be small numerically, we introduce the ratio y* 
that measures the relative change of y* with increasing 
Eb, 


where yg is the value of y* for Eb = 0 (i.e., the fully flex¬ 
ible case). Apparently, a positive or negative y* implies 
that chain stiffness promotes or opposes self-assembly, 
respectively. 

Figure 2 displays y* as a function of Eb for various 
sticky interaction energies eg, when the other parameters 
are held constant, y* hrst increases with Eb and then sat¬ 
urates at sufficiently large Eb for small |es|. Moreover, the 
quantitative dependence on Eb weakens as |es| increases 
and even disappears when |es| is sufficiently large (e.g., 
|es| = 3000 K), and the saturated value for y* decreases 
with increasing |es|. This is understandable because y* 
is shown to increase with |es| and reach its upper limit 
(i.e., 4>/M) for large Clearly, y* fails to elevate 

with increasing Eb after y* attains its upper limit. Fig¬ 
ure 3 exhibits y* as a function of Eb for various molecular 
weights M. As can be seen, y* grows more rapidly for 
larger M and the saturated value for y* increases with 
M. For example, the saturated limit for y* grows from 
0.023 to 0.079 when M is increased from 5 to 100. Of 
course, the saturated value for y* remains independent of 
M after M reaches a large value above which most poly¬ 
mer properties become insensitive to changes in molec¬ 
ular weight. Therefore, Figs. 2 and 3 indicate that the 
influence of chain stiffness becomes more significant for 
telechelic polymers with weaker sticky interactions and 
higher molecular weights. 

While Figs. 2 and 3 may suggest that the chain stiff¬ 
ness generally promotes self-assembly. Fig. 4 indicates 
that the model of linear telcehelic polymer melts per¬ 
mits self-assembly also to be suppressed by chain stiff¬ 
ness, a trend that occurs when (j) is sufficiently low at 
T = 200 K. Further analysis indicates that temperature 
also bears on whether chain stiffness promotes or opposes 
the self-assembly. For example. Fig. 5 displays the con¬ 
tour plot of y* in the 4>-T plane, where the bending en¬ 
ergy is fixed to be Eb = 2000 K. Although chain stiffness 
promotes self-assembly for systems represented in the (fh 
T plane where either ^ or T is high (termed region I), 
self-assembly is indeed suppressed by chain stiffness when 
both (j) and T are sufficiently low (termed region II). 

The result in Fig. 5 may be rationalized as follows. 
Linear clusters are expected to form more easily than 
cyclic clusters (i.e., rings) at high volume fractions, in 
line with a FH type theory of self-assemblyi^ The pic¬ 
ture clearly changes when polymer volume fraction is low. 
In particular, the FH type theory of self-assembly2£ pre¬ 
dicts that rings predominate over the formed clusters at 
low temperatures, whereas the opposite situation ensues 
at high temperatures, a behavior that arises because the 
extra bond energy gained upon ring closure outweighs the 
entropy loss upon ring closure as temperature decreases. 
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FIG. 6. Dependence of the boundary line in the (f)-T plane on (a) the sticky interaction energy ts, (b) the bending energy Eb, 
(c) the molecular weight M, and (d) the van der Waals interaction energy e. 


Therefore, the formation of cyclic clusters is favored only 
when both <f) and T are low; otherwise, most of the clus¬ 
ters that form in the system are linear. The trend of 
forming linear clusters is enhanced by chain stiffness be¬ 
cause of a diminished probability of ring closure as the 
chains stiffen. Moreover, the formation of sticky bonds 
between different chains becomes less sterically hindered 
as the chains stiffen. Hence, it is observed in Fig. 5 
that chain stiffness promotes self-assembly in region I, 
where linear clusters predominate. Meanwhile, if cyclic 
clusters predominate, the gain in sticky bonds due to the 
enhancement of linear clusters induced by chain stiffness 
cannot compensate for the loss of sticky bonds due to 
the stiffness generated reduction in ring formation. Con¬ 
sequently, chain stiffness opposes the self-assembly, re¬ 
sulting in the appearance of region II in the (/)-T plane in 
Fig. 5. Clearly, the above explanation is largely based 
on a FH type theory;^ because the present LCT pro¬ 
vides no information concerning the formation of rings. 
On the other hand, computer simulations may provide a 
clearer microscopic picture for understanding the trends 
for the variation of the average degree of self-assembly 
with chain stiffness. 

The presence in Fig. 5 of two regions with opposite 
behavior implies the existence of a boundary in the (jyT 
plane that has y* = 0 and, therefore, that separates the 
two regions with opposite dependences of y*on chain stiff¬ 


ness (see the solid line in Fig. 5). Of course, the concept 
of a boundary line in the (p-T plane becomes meaning¬ 
less for linear telechelic polymers with sufficiently large 
|es|, where chain stiffness ceases to affect y* (see Fig. 2). 
An examination of the variation of the boundary with 
molecular parameters reveals that the boundary line is 
insensitive to Cg [Fig. 6(a)]. The area of region II shrinks 
slightly with increasing Ei, and saturates for sufficiently 
large E}, [Fig. 6(b)], a trend that ensues upon increasing 
M, as shown in Fig. 6(c). Moreover, Figure 6(d) indi¬ 
cates that elevating e leads to a dramatic increase in the 
area of region II in the 4>-T plane. 


B. Influence of chain stiffness on the enthalpy and 
entropy of self-assembly 

The enthalpy Ahp and entropy Asp of self-assembly 
characterize the thermodynamics of self-assembly in clas¬ 
sical FH type theories for the equilibrium self-assembly 
of linear polymers^^^— but they are commonly treated 
as phenomenological parameters. On the other hand, the 
LCT for telechelic polymers allows the determination of 
Ahp and Asp as a function of thermodynamic and molec¬ 
ular parameters. The previous analysis^^ indicates that 
Ahp and Asp strongly depend on various thermodynamic 
and molecular parameters, and a similar analysis shows 
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FIG. 7. Enthalpy Ahp and entropy Asp of self-assembly of 
linear telechelic melts as a function of bending energy Ei, for 
various polymer volume fractions 4>- The computations are 
performed at a constant temperature of T = 300 K for melts 
of linear telechelic chains, where the molecular weight is M = 
5, the van der Waals interaction energy is e = 100 K, and the 
sticky interaction energy is ts = —300 K. 


tionsr^ 

/3A/ip = /3^ 

and 


1 d[Pfs) 


^y* 9/3 




= /3e.-^(y*ri/3^(7) 


ks y* I ks dT 
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ks 
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where Sg^comb is the portion of the system’s total combi¬ 
natorial entropy that depends on the density y* of sticky 
bonds, 


ks 


=(j)X ln((/)a:) — {(j)x — 2y*) h\{(j)x — 2y*) 


-y 


l + lnlHf 


(9) 


In the limit of infinite dilution (i.e., ^ = 0 and (j) = l)i 
Eqs. (7) and (8) simplify significantly for fully flexible 
chains,— 

Ahp{Eb 0) = Cs - (1 -b -)e,for (j) = 0, (10) 


Ahp{Eb ->• 0) = es,for (j) = 1, 


( 11 ) 


Asp{Eb —>■ 0) 2 5 


ks 


= -2,for(/) = 0, (12) 




that chain stiffness can be added to the list of molec¬ 
ular properties that greatly influence the enthalpy and 
entropy of self-assembly. 

Because the present theory employs a coarse grained 
model in which each united atom group effectively com¬ 
prises, say, an actual monomer, changes in enthalpies 
and entropies on association in real systems may have 
contributions from the missing internal degrees of free¬ 
dom that would have to be included as parameters in the 
present formulation but that would be derived by includ¬ 
ing these additional degrees of freedom by describing the 
telechelics with structured monomers. When studying 
systems with solvent molecules, there may be additional 
contributions from solvent reorganization, disruption of 
solvent network, etc. The present theory neglects the 
contributions arising from the internal degrees of free¬ 
dom. 
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and for fully stiff chains, 

Ahp(Eb —>• oo) = Cs — e, for <() = 0, 


(13) 


(14) 


Ahp{Eb —>• oo) = Cs, for <() = 1, 


(15) 


The enthalpy Ahp and entropy Asp of self-assembly 
are defined as the derivatives of the portion fg of the 
free energy / that arises entirely from the sticky interac- 
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(17) 


As shown in Eqs. (10-17), the enthalpy Ahp and en¬ 
tropy Asp of self-assembly display considerable variation 
even in the limit of infinite dilution when the chains 
shift from fully flexible to fully stiff. For example, Ahp 
and Asp increase by {2/z)e and 6/z'^ in the limit of 
(p = 0, respectively, when the bending energy elevates 
from Eb 0 to Eb ^ oo. While both limits {Eb —>■ 0 
and Eb —?> oo) lead to an identical value of Ahp = for 
(j) = 1, a finite Eb leads to a decrease in Ahp ioi p = 1 
[see Fig. 7(a)]. The entropy Asp of self-assembly in the 
limit oi Eb ^ oo is generally found to be larger than that 
in the limit of Fb —0. Figure 7 also implies that the de¬ 
pendence of Ahp and Asp on Eb changes sensitively with 
other parameters such as ip. 


C. Influence of chain stiffness on the entropy density of 
linear telechelic melts 


As expected, the chain stiffness of linear telechelic 
polymer melts exerts a major impact on the entropy den¬ 
sity (i.e., the entropy per lattice site), which is a central 
quantity in the GET of polymer glass formationJ^ The 
analytic expression for the entropy density simply follows 
from the specific Helmholtz free energy / as 


s = — 



(18) 


Likewise, the entropy density is the sum of two separated 
terms So and Ss , corresponding to the contributions aris¬ 
ing from the reference system and the sticky interactions, 
respectively. Consequently, So and Ss are defined by 


So = - 


dfo 

dT 


(19) 


and 



( 20 ) 


Following previous work^i^ the computations are per¬ 
formed at constant pressure P, which is determined from 
the Helmholtz free energy F, 


P = - 


dF 


1 OF 


m,T 


Fcell dNi 


( 21 ) 
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FIG. 8. Entropy density s/fcs and its various contributions 
SojkB and Ss/ks as a function of temperature T for a melt 
of fully flexible linear telechelic chains (red lines) and a melt 
of semiflexible linear telechelic chains (green lines) for (a) a 
weak sticky interaction strength ts = —300 K and (b) a strong 
sticky interaction strength = —3000 K. The computations 
are performed at a constant pressure of P = 1 atm. The 
molecular weight is M = 8, and the van der Waals interaction 
energy is e = 200 K for both melts, while the bending energy 
is Eb = 400 K for the semiflexible melt. The dotted line in 
(b) highlights the fact that the contribution from the sticky 
interactions to the entropy density can be negative for very 
strong sticky interactions. 


where V is the volume of the system, and Well = Oceii is 
the volume associated with a single lattice site. The pa¬ 
rameters z = 6, Oceii = 2 . 7 A, and e = 200 K are used in 
our illustrative computations. A low molecular weight of 
M = 8 is chosen for the examples because the contribu¬ 
tion from the sticky interactions to the entropy density 
becomes more significant for lower M ; see discussion be¬ 
low. 

Figure 8 displays the temperature dependence of the 
entropy density s as well as its various contributions So 
and Sg for a melt of semiflexible linear telechelic chains 
with the bending energy Eb = 400 K at a constant pres¬ 
sure of P = 1 atm for weak sticky interactions = —300 
K [Fig. 8(a)] and strong sticky interactions Cg = —3000 
K [Fig. 8(b)], respectively. The variation for fully flex- 




























ible linear telechelic polymer melts is also included for 
comparison, as well as to demonstrate the importance of 
including chain semiflexibilty within the LCT. (The ex¬ 
pression for the specific free energy for fully flexible linear 
telechelic polymer melts can be found in Ref. [ll|.) The 
term So from the reference system dominates the entropy 
density s of the system for both weak and strong sticky 
interactions. In fact, the contributions arising from the 
sticky interactions become even less significant for larger 
molecular weights, and s is nearly identical to Sq for suf¬ 
ficiently large M. This behavior agrees with expectation 
since the concentration y* of sticky bonds becomes ex¬ 
tremely small for very large molecular weights Recall 
that the upper limit for y* is inversely proportional to M 
as y^ax = 4>/M. The above feature applies for both fully 
flexible and semiflexible linear telechelic polymer melts. 

Furthermore, we find that increasing the bending en¬ 
ergy has no significant influence on Ss but leads to a 
dramatic drop in So (see either panel of Fig. 8). Conse¬ 
quently, the entropy density s for a semiflexible melt de¬ 
creases precipitously upon cooling and approaches zero at 
low temperatures, leading to the appearance of a charac¬ 
teristic feature of glass formation within the GET. (The 
vanishing of the entropy density is probably an artifact 
of the high T expansion in the LCTi^) Therefore, it is 
possible to explore glass formation in telechelic polymers 
using the GET only if the polymer chains are modeled as 
being semiflexible, although the fully flexible telechelic 
model suffices in establishing the trends obeyed by the 
thermodynamic properties that are weakly dependent on 
chain stiffnessi ^^d^ (The maximum displayed by the en¬ 
tropy density s in Fig. 8 is also a characteristic feature 
of glass formation within the GET, and Ref, [isl discusses 
why the entropy density depends non-monotonically on 
temperature.) Eigure 8 also reveals that the term aris¬ 
ing from the sticky interactions can be either positive or 
negative, depending on whether the sticky interactions 
are weak or strong. Thus, it suggests that the sticky 
interaction energy Cg may have interesting influences on 
glass formation in linear telechelic polymer melts. Fu¬ 
ture work will combine the current extension of the LGT 
with the AG relatiou j^^i^^ to provide a generalization of 
the GET that enables detailed analysis for the role of 
self-assembly and the magnitude of the sticky interaction 
energy Cg upon glass formation in telechelic polymers. 


IV. SUMMARY 

We apply the LGT for semiflexible linear telechelic 
polymers to assess the influence of chain stiffness on 
the basic thermodynamic properties of self-assembling 
telechelic polymer melts and specifically illustrate general 
trends for the dependence of the average degree of self- 
assembly upon chain stiffness. The calculations from the 
LGT imply that chain stiffness promotes self-assembly of 
linear telechelic polymer melts when either the polymer 
volume fraction or the temperature is high, but opposes 


the self-assembly when both polymer volume fraction and 
temperature are sufficiently low. A FH type theory^ of 
the competition between the formation of rings versus lin¬ 
ear clusters is invoked to provide a possible rationale for 
the predictions from the LGT concerning the influence of 
chain stiffness on the average degree of self-assembly. Our 
results indicate that chain stiffness provides an important 
molecular variable for tailoring the physical properties of 
self-assembling telechelic polymers. Meanwhile, we em¬ 
phasize that simulations clearly offer the possibility of 
better understanding the influence of chain stiffness on 
the average degree of self-assembly in telechelic polymers 
by analyzing the structural and dynamical properties, 
which are inaccessible by our thermodynamic model. Our 
theory instead provides a tool for guiding the design of 
telechelic polymer materials by establishing the relation 
between the molecular details and the thermodynamic 
properties. 


FH type theories for the equilibrium self-assembly 
of polymers traditionally employ highly coarse grained 
models)^^— where the interaction parameters (such as 
the enthalpy and entropy of self-assembly) that char¬ 
acterize the thermodynamics of self-assembly must be 
adjusted phenomenologically. The LGT for telechelic 
polymers instead provides a theoretical tool for deter¬ 
mining these interaction parameters as a function of 
molecular and thermodynamic parameters of the self¬ 
assembling system. While a previous paper— illustrates 
the strong dependence of the enthalpy and entropy of 
self-assembly on temperature, polymer volume fraction, 
molecular weight, van der Waals interaction energy, and 
sticky interaction energy, the present paper continues to 
demonstrate the significant influence of the bending en¬ 
ergy on these interaction variables. 


One substantial benefit of the LGT for semiflexible lin¬ 
ear telechelic polymers lies in the fact that the charac¬ 
teristic features of glass formation appear in the GET— 
only if the polymer chains are modeled as being semi¬ 
flexible. Hence, we provide illustrative calculations for 
the influence of chain stiffness on the entropy density 
of linear telechelic polymer melts and demonstrate the 
importance of including semiflexibility within the LGT 
for exploring glass formation in linear telechelic polymer 
melts. Our illustrative calculations also imply that the 
contribution arising from the sticky interactions to the 
entropy density can be either positive or negative, de¬ 
pending on whether the sticky interactions are weak or 
strong, suggesting that the sticky interaction energy may 
have interesting influences on glass formation in linear 
telechelic polymer melts. A generalization of the GET 
for exploring the influence of self-assembly on glass for¬ 
mation in linear telechelic polymer melts can be achieved 
by combining the current extension of the LGT with the 
AG relation— 
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